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Abstract: We study the supersymmetric properties of D7-branes in the warped deformed 
conifold. We consider the K-symmetry conditions on D7-branes in this specific warped back- 
ground, taking into account the background NS-NS 2-form flux. While any holomorphic 
embedding defines a supersymmetric D7-brane in the absence of background H-flux, most 
i-£^ of the D7-brane embeddings considered in the literature do not preserve supersymmetry 

for the warped deformed conifold without also including brane worldvolume flux. For the 
simplest such embedding, we construct numerically the worldvolume flux necessary to re- 
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store supersymmetry. We also comment on the dual field theory descriptions in terms of 
^ l cascading J\f = 1 supersymmetric gauge theories with flavors. Finally, we discuss some 

C\J possible applications of our results to moduli stabilization and vacuum energy uplifting, 

f — gauge/gravity duality, and string inflationary model building. 
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1. Introduction 

One of the major developments in string theory of the last decade has been the emergence 
of warped spacetimes as a setting for a rich variety of physics. To date, perhaps the most 
important application of such warped geometries has been in the gauge theory/string theory 
correspondence [1, 2, 3]. In this setting, noncompact spacetimes which are warped by the 



presence of background flux have been conjectured to give a dual description of gauge 
theories. In many cases, this duality can relate strongly coupled field theories to weakly 
coupled gravity theories, and as such offers new tools for studying field theories when 
traditional perturbative methods fail. 

Another application of warped geometry has been in the study of string theory as a 
theory of particle physics. When one compactifies string theory from ten dimensions to 
four dimensions, the compactification manifold can support nontrivial background fluxes, 
and these fluxes will generically have a backreaction on the geometry, resulting in warping. 
The gravitational redshift generated by a strong warp factor can lead to a hierarchy of 
scales, thus offers a geometrical explanation for the huge disparity between the electroweak 
and the Planck scale [4] (see also [5, 6, 7, 8] for string theory realizations). For the same 
reason, warped geometries have also been invoked in supersymmetry breaking scenarios 
in string theory, both to lower the scale of supersymmetry breaking [9, 10] and to control 
the amount of sequestering [11]. The low energy effective supergravity theory describing 
strongly warped backgrounds is challenging to derive, though continued progress has been 
made [9, 10, 12, 13, 14, 15, 16]. 

Warping has also found application in inflationary cosmology. A key issue in infla- 
tionary cosmology is to understand the underlying dynamics which allows for a sufficient 
number e-folds of expansion. It turns out that warping can be fruitfully used in stringy 
cosmological models to construct suitable potentials for slow roll inflation as well as for 
motivating new inflationary mechanisms (see [17] and the references therein). 

Strikingly, there is a particular warped geometry, the warped deformed conifold in 
type IIB supergravity [18, 19], which is simple enough to allow detailed study, but is 
also rich enough to have been used in all three of these areas. It has almost anti-de Sitter 
asymptotics, and in gauge/gravity duality it is dual to an interesting confining gauge theory. 
If one embeds the deformed conifold in a compact Calabi-Yau, and turns on appropriate 
three-form fluxes, one finds that the resulting background has a natural and stable hierarchy 
of scales [8] . And in cosmology, interesting inflationary models have been also constructed 
in the conifold [20, 21, 22, 23]. 

In this paper, we study a particularly interesting variation of the warped deformed 
conifold, in which we add a certain number of probe D7-branes to the background. These 
branes fill four noncompact directions and wrap a four-cycle £4 in the transverse space. 
In the gauge theory dual, these branes add fundamental matter to the theory [24, 25], 
offering the prospect of field theory duals which resemble QCD. Moreover, in more phe- 
nomenological applications, the D7-branes allow the introduction of gaugino condensates, 
which can stabilize Kahler moduli [26] and which can introduce potentials for mobile D3- 
branes (candidate inflatons in the scenarios considered in [20, 21, 22]). Given this range of 
applications, it is of clear interest to find as many supersymmetric D7-brane embeddings 
in the deformed conifold as possible, allowing for greater freedom in model-building. In 
this paper we study the supersymmetry conditions for several different classes of branes, 
focusing in particular on the cases where worldvolume flux is necessary for the brane to 
preserve supersymmetry. 

This paper is organized as follows. We begin by reviewing the deformed conifold in 



Section 2, along the way collecting some results which will be useful in the rest of the 
paper. In Section 3, we review the criteria for the D7-branes to preserve supersymmetry 
(the K-symmetry conditions) in a general Calabi-Yau flux compactification. We make 
use of these results in Section 4, where we study the case of D7-branes in the conifold 
with background fluxes. There is one such D7-brane embedding (sometimes known as the 
Kuperstein embedding) already known to be supersymmetric with no worldvolume flux 
turned on [27] . On the other hand, the embedding first studied in [28] is not supersymmetric 
when the worldvolume flux vanishes. Our main result is in Section 4.3, where we construct 
the worldvolume flux necessary to restore supersymmetry, by a combination of analytic 
and numerical methods. We also offer a criterion that in many cases can show whether a 
worldvolume flux is needed to restore supersymmetry, and which is far simpler to verify 
than explicitly checking /{-symmetry. In Section 5, we make some observations about the 
field theories dual to the conifold with our D7-branes included. Finally, we end with a 
discussion of prospects for the future in Section 6. 

2. Review of the Warped Deformed Conifold 

In this paper we focus on a particular concrete example of a warped Calabi-Yau space, 
namely the warped deformed conifold. This geometry is noncompact, but in principle it 
can be embedded in a compact Calabi-Yau with background bulk fluxes. We now proceed 
to review some facts about the conifold and warped compactifications in general which will 
be useful in the rest of the paper. 

The conifold is most simply defined as a submanifold of flat C 4 given by the equation 
[18] 

z\ + zl + zl + z\ = Q (2.1) 

where the Z{ are coordinates of the ambient C 4 . This space has a manifest £77(2) x 
577(2) x U(l) global symmetry, where the 5(7(2) x 577(2) ~ SO (4) acts by rotating the Z{ 
and the U(1)r acts as multiplication of each of the Zi by a phase, Z{ — > e tip Zi. The defining 
equation of the conifold is also invariant under an overall rescaling of the coordinates, 
which implies that it admits a conical metric. The corresponding geometry is a cone over 
a five-dimensional Einstein manifold called T ' and the metric takes the form 

dsl = dr 2 + r 2 ds^ 1}1 . (2.2) 

where the metric of the T 1 ' 1 base is 

d4u = g ( # + Yl cos 6id & ) + 6 ^ ^ + sin2 ° id $) (2 - 3) 

and the angular coordinates range as < ip < 4ir, < 6% < ir, and < 0, < 2ir. From the 
form of the metric it is clear that base is a sphere fibration S 1 x S 2 x S 2 and it turns out 
that its topology is actually S 3 x S 2 . The base space has Betti numbers 62 = ^3 = 1, and 



the associated harmonic forms are 

w 2 = \ (On - 22 ) (2.4) 

W3 = (Aw 2 (2.5) 

where £ = dtp + cos d\d§\ + cos #2^2 an d 

fijj = d^j A sin 6 jd(j)j. (2-6) 

This space is singular at the tip of the cone where all the 2$ are zero (and so we will 
sometimes refer to this space as the "singular conifold.") 

It is also possible to define the conifold in terms of another set of complex coordinates, 
by the equation 

W1W2 — W3W4 = (2-7) 

which is related to the original defining equation by an obvious change of variables. It 
turns out that this alternative description will also be useful to us. One nice property is 
that the W{ coordinates admit a relatively simple description in terms of the angles on T > : 

W1 = ^w-fc-ao sin 0i sin 02 (28) 

W2 = f 2/* e */W+*i+<h) cos °1 cos |, (2.9) 

W3 = r 3/2 e W-* 1+ *0 sin °1 cos fc (2J0) 

W4 = r 3/2 e ,/2(^ 1 -^) cog °1 gin °1, (2 .11) 

One of the important properties of the conifold is that its singularity can be smoothed 
by adjusting moduli. At the tip of the cone, one may perform either a small resolution, 
blowing up an S 2 , or a deformation, by changing the defining equation of the conifold to 

z\ + zl + zl + zl = e 2 . (2.12) 

We can make a phase rotation of the z% so that e is real, and we will always assume this 
in the rest of the paper. It is evident from taking a real slice of this defining equation 
that deformation results in the appearance of a finite-sized S s at the tip of the cone. The 
deformed conifold preserves the <SO(4) isometry of the singular conifold, but the scale 
invariance is broken by e and the U(1)r symmetry is broken to a Z2 subgroup. 

This deformed geometry has special relevance in flux compactifications because of its 
topology. The existence of a three-cycle allows us to turn on RR 3-form flux on this cycle, 
and there are supersymmetric supergravity solutions if we also include NS-NS flux on the 
(noncompact) dual cycle. 

We now proceed to review flux compactification of the type considered in [8] , where 
a warped deformed conifold throat naturally develops near the local conifold singularities. 



Let us begin with the ten dimensional warped metric preserving four-dimensional Poincare 
symmetry which takes the following form 

ds 2 1Q = e 2A ^~g^dx»dx» + e- 2A ^g mn dy m dy n . (2.13) 

Here g^ ,fj,,i> = 0,1,2,3 and g mn ,m,n = 4,... ,9 are the unwarped metrics for the 
non-compact four dimensional space time and the compact six dimensional manifold A4q 
respectively, whereas e ( y > is the warp factor which only varies over Mq. 

Following [8] , we shall allow for both RR and NS-NS fluxes of IIB supergravity to be 
turned on, and they can be written succinctly as: 

G 3 = F 3 - tH 3 = \G mnp [dy m A dy n A dy p ] , (2.14) 

D 

F 5 = F 5 - -C 2 A H 3 + -B 2 A F 3 = (1 + *i Q )[da(y) A dx° A dx 1 A dx 2 A dx 3 } . (2.15) 

Here we have allowed for IIB axio-dilaton r = r{y) = Co — ie~^ to vary over M§ and 
in (2.15) we have used self-duality F§ = *\qF§, four dimensional Poincare invariance and 
Bianchi identity to constrain the expression for F§. In addition, local objects extending 
into non-compact four dimensions can also be included (which may wrap cycles in Me), 
and they need to satisfy the tadpole cancellation condition: 

-4- f H 3 AF 3 + Q l £^ = 0, (2.16) 

where Q£ is the D3-brane charge on the local objects. The supergravity equation of 
motion for such configuration of fluxes and local sources then yields: 
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(2.17) 



where we have defined the following linear combination: 

$± = e 4A ^ ± a(y) , G± = iG 3 ± * 6 G 3 , (2.18) 

and T^ u and T mn are four and six dimensional energy-momentum tensors. A special solution 
to (2.17) is given by the so-called imaginary self-dual (ISD) configuration: 

$_=0, G_=0, <-► e 4A{y) = a{y) , * 6 G 3 = iG 3 , (2.19) 

(rpm t^M\ local 

1 m /*> = T 3/ 4 ocal . (2.20) 

In terms of the local sources in IIB, the equality (2.20) can be satisfied by D3 branes or 
03 planes, D5 brane wraps on vanishing two cycle, also known as "fractional D3"; or D7 
wraps on non-trivial four cycles, e.g. K3. Turning now from these generalities back to the 
conifold, let us turn on M units of F 3 along the so-called A-cycle which corresponds to the 
minimal S 3 in the deformed conifold, and also (-K) units of H 3 in the dual B-cycle. The 
quantization condition for the fluxes then yields: 

— — [ F 3 = M, r / H 3 = -K. (2.21) 

2ira> J A ' 2ira> Jb 



These three form fluxes stabilize the value of e 2 at 

e 2 ~ exp (-2vr K/g s M) , (2.22) 

and the warp factor at the tip is given by e ° = e - 2 ^ K / 3M 9s The total D3 charge is given 
by the intersection form: 

H 3 AF 3 = MK = N . (2.23) 
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The explicit IIB supergravity solution corresponds to the flux configuration described here 
was given in [19], where the authors started with N unit of D3 branes and M unit of 
fractional D3 branes at the tip of singular conifold, with g s N,g s M S> 1. In the near 
horizon limit, these branes are replaced by their fluxes, one can obtain explicit expressions 
for the corresponding warp factor and fluxes (for reviews, see e.g. [29], [30]). 

In the large radius limit, the supergravity solution including these 3- form fluxes (but 
with the axion and dilaton constant) has the useful description (due to Klebanov and 
Tseytlin, or KT) [31]: 

ds\ Q = e 2A ^(d Xfl dx^) + e~ 2A ^ {dr 2 + ^ds 2 ^) (2.24) 

~ AA{V) = §s*9^ (n + ^-g s M 2 log(r/r )) (2.25) 

9s F 5 = d 4 x A de 4A(r) + * w (d*x A cie 4A(r) ) (2.26) 

Ma' 
F 3 = —lu 3 (2.27) 

B 2 = -^ log(r/r ) A u 2 (2.28) 

„ 3o,Afa' dr 
H 3 = dB 2 = -^- A W2. 2.29 

2 r 

This solution possesses a naked singularity at small r; however, in the full KS solution the 
deformation of the conifold resolves the singularity. 

For later use, we will record here the explicit Kahler form and the NS-NS B-field on 
the warped deformed conifold. 50(4) invariant expressions can be found in [29]: 

J = K'{p) l^2dz z A dzi j + K"(p) I J2 ^ dz i ) A I 2 z i d *i ' ( 2 - 30 ) 

„ , , s _ , ,_ ,, s ig s Ma' pcothp — 1 , n n -,\ 

B 2 = b{p)e ijk iz i z j dz k A dz x , b{p) = — -* , (2.31) 

2|ep sinn p 

where KL{p) the Kahler potential for the deformed conifold is defined implicitly as: 

f sinh(2p) - 2p \ 5 

K {p) = VW^7) ' ( ] 

K,"(p) = ( 5 + cosh(2p) - 6p coth p) ^ 
6 1 e | 2 sinh p(p — cosh p sinh p) 



Here () denotes the differentiation with respect to r 3 , and the usual radial coordinate r is 
related to the coordinate p by 



r 3 



^\ Zl \ 2 = |e| 2 coshp. (2.34) 
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The following definitions, taken from [30], prove to be notationally convenient: 

4 4 4 

m = e ijk iZiZjdz k A dzi , 7? 4 = 'Y^z i dz i A^Zjdzj , rj 5 = ^ dzi A dz t , (2.35) 

i=l j=l i=l 

and we will make use of these forms in the calculations to follow. In the large p limit, b{p) ~ 
pe -4 exp(— 2p) ~ r _6 log(r) and one finds that the potential Bi varies logarithmically with 
r, as in [31]. 

3. Supersymmetry Conditions 

Given the rich physics that has arisen from the study of D7-branes in warped flux 
compactifications, it is interesting to find as many different stable configurations of D7- 
branes in these backgrounds as is possible. Because supersymmetry guarantees stability, 
and is typically a simpler property to check than non-supersymmetric stability, we will 
focus on a search for supersymmetric D7-branes in warped throats with flux. We begin 
by reviewing the supersymmetric embedding conditions given in [32, 33] for D7 branes 
wrapping a four-cycle in a generic Calabi-Yau three-fold, and then in the next section we 
will apply these conditions to specific examples of embeddings in the warped deformed 
conifold. 

Consider a CY three-fold A4q with Kahler two-form J, a background NS-NS 2-form 
potential B 2 , and some number of coincident D7-branes which fill out the four transverse 
Minkowski space-time dimensions and extend in four directions on the compact manifold. 
There is a vector field A supported on the D7-brane worldvolume, and its curvature F% = dA 
combines with the background B 2 to form a gauge invariant combination J-. 

F = B 2 + 2WF2 , (3.1) 

where we use " to denote the pull-back of bulk quantities on to the D7-brane worldvolume. 
In [32, 33], it was found that the K-symmetry condition for D7 branes wrapping a 
four-cycle S4 in Ai§ are as follows: 

1. The four cycle S4 must be holomorphic. In other words, the D7-brane locus may be 
written as a holomorphic equation in the complex structure of the CY. 

2. The generalized J- field strength is a two-form of pure (1, 1) type, or 

jC-2,0 = jc-0,2 = . (3.2) 



3. The Kahler form J pulled back to S4 and the generalized field strength T satisfy the 

equation: 

a ~ fe 2A ~ 1 - A 

e- A JAF = tan9(—JAJ--FAF). (3.3) 

In (3.3) we have included the warp factor e in (2.13) and the constant 9 is specified by 
the supersymmetry conditions for a given supergravity background solution [32]. These 
conditions are equivalent to the requirement that the D7-brane configuration is a local 
energy minimum, and in the fluxless case are equivalent to the condition that E4 is a 
minimal surface. 

J- must also satisfy a modified Bianchi identity: 

df = H 3 , (3.4) 

where H3 is the pull-back of H3 = dE>2- Equation (3.4) arises from the definition of j- 
when the gauge field strength by itself is closed and is necessary to guarantee cancellation 
of the global anomaly for open fundamental strings ending on the D7 branes [34]. Taken 
together, the condition (3.4) and the supersymmetry conditions are sufficient conditions 
for the equations of motion from the D7 brane action to be satisfied. 

In the case of the warped deformed conifold, the parameter 6 is zero and the condition 
(3.3) simplifies to 

JAf = 0, (3.5) 

JF satisfying (3.5) is said to be "primitive." An alternative and very useful interpretation 
[32] for conditions (3.2) and (3.5) is that the generalized field strength JF is anti-selfdual: 

jr=-* 4 jr ) (3. 6 ) 

where the Hodge star operation *4 is taken with respect to the induced metric on the four 
cycle. 

In terms of complex embedding coordinates {z{}, J and B2 are manifestly (1,1). There- 
fore in the absence of the magnetic field F2, T is also of the type (1,1) satisfying the 
supersymmetry condition Eqn. (3.2). 

4. Supersymmetric D7 Embeddings in the Conifold 

In this section we shall examine several holomorphic D7-brane embeddings in the singular 
conifold that have appeared in the literature. In the absence of background H-flux, the 
only requirement for a D 7- brane to be supersymmetric is holomorphy of its embedding 
equation. Several such holomorphic embeddings have been proposed, and in the singular 
conifold limit, the K-symmetry conditions are trivially satisfied, and in some cases the 
supercharges have been explicitly constructed [37]. However, the supersymmetry conditions 
are more complicated when H-flux is turned on, and our goal in this section is to explore 
this situation in detail. 



We will focus first on linear holomorphic embeddings in the interest of simplicity. 
There are two natural classes of linear embeddings, given by either the description of the 
conifold in terms of the 2j coordinates or the Wi coordinates. The embedding z\ = fi was 
studied in detail in [27] and it satisfies the supersymmetry conditions (3.5) even on the 
deformed conifold without needing any additional world volume flux. Curiously, among 
the commonly studied classes of embeddings it appears to be the only one that does not 
require a flux. 

Another simple class of linear embeddings was introduced in [28] and is most naturally 
expressed in the Wi coordinates as, for example, w\ = //. When nontrivial H-flux is turned 
on, this embedding is not supersymmetric without also turning on the world volume field 
strength F 2 . We will exhibit the necessary flux, thus demonstrating that this class of 
D7-branes can in fact be made supersymmetric. To make the calculations tractable we 
have taken the large-radius (KT) limit of the warped deformed conifold, but it would be 
interesting to study these branes in the full KS geometry as well. 

In the following section we will refer to the relevant D7-branes as being in either 
the "2-embedding" or the "u>-embedding" , depending on the natural coordinates of the 
embedding equation. Having understood the simplest D7-brane embeddings, we will also 
see that many other more complicated embeddings will require worldvolume flux to be 
supersymmetric . 

4.1 2-Embedding 

We begin by considering the embedding for the deformed conifold proposed in [27] 

Zi = fi, /ieC. (4.1) 

This embedding breaks the SO (4) isometry group of the deformed conifold down to the 
50(3) which rotates {21,22,23}. It was shown to be supersymmetric in the deformed 
conifold for real // [27]. We review the calculation here, generalizing as well to complex //. 

In order to check the supersymmetric condition (3.3), the pull back of the NS-NS B-field 
for this embedding can be obtained from 

Vi = »7i|*4=a* = Biidzi A dz\ + B l2 dz\ A dz 2 + B 2 \dz 2 A dz~\ + B 22 dz 2 A dz 2 , (4.2) 

where the various components are given by: 

B H = -\Z3\~ 2 (ZIZ 3 - Z 3 Z 1 ){z 2 Jl - Z 2 fJ.) , 

B\2 = ( z 3fi ~ zz^) + |z 3 |" 2 (zizz{zip, - z\y) + z 2 z 3 (z 2 fi - Z 2 fi)) , 

B 2 \ = -(Z3ft - Zzfj) - \Z3\~ 2 (ziZ 3 (zifl - 21//) + Z 2 Z 3 (z 2 fl - 22//)) , 

B 22 = \z 3 \~ 2 (z 2 z 3 - z 3 z 2 )(z 1 Jl - 21//) , 



and 



23 = ^(e2 - ^j - (zf + zl) . (4.3) 



The components of the Kahler two-forms J are given by: 

V5 = ifo\(zi=n) = J[\ dz\ A dz\ + J' X2 - dz\ A dz 2 + J' 2 i dz 2 A dz\ + J^ dz 2 A dz 2 , (4.4) 
774 = r/4|( Z4=M ) = J"j dz\ A dzi + J"g dz\ A dz2 + «/ 2 I dzi A etei + J 22 - dz 2 A <iz2 • (4-5) 

Here the various functions are given by: 

ji _ 1 , N 2 7/ _ 1 1 N 2 ./ _ ^1^2 r _ Z1Z2 ,. n 

11 + U Q |2 ' J 22 _ i_h | r l2> J 12 _ U„|2 ' J 21 _ | |2 - V 4 -°J 

|^o I |^3| |^3| |*d I 

T // __ (jjjjj - ^3^l) 2 T „ __ (z 2 ^3 ~ Z3Z2) 2 T „ _ -,/ __ (Z1Z3 ~ Zj,Zi){z 2 Zj, ~ Z3Z2) 

11 ~ UJ2 ' J 22 - I 12 ' J 12 - J 21 - I 12 

1^31 I ^3 1 |^3| 

(4.7) 
By either tedious or computer assisted algebra, one now obtains 

V5 A 771 = , (4.8) 

774 A 771 = . (4.9) 

From equations (4.8) and (4.9), we can confirm that the SUSY condition (3.3) is satisfied 
by the holomorphic embedding (4.1), as noticed in [27]. Moreover as T = B 2 , the Bianchi 
identity (3.4) is also trivially satisfied. 

4.2 ^-Embedding 

Let us now consider the holomorphic embedding given in [28]: 

w 1 = Z ^^= f i, /iGC. (4.10) 

v2 

This embedding breaks the 50(4) = SU(2) x SU(2) rotational symmetry of the deformed 
conifold to the U(l) subgroup which rotates w 3 and w± by opposite phases. By combining 
(4.10) with the defining equation of the deformed conifold, we can also deduce that 

w = ^ 5 = ^M ±ii ) i (4u) 

which leaves us only {z 3 , 214} or {w 3 , W4} and their complex conjugates as the independent 
variables parameterizing the four cycle wrapped by the D7 brane. In this case we shall first 
demonstrate that supersymmetry is broken and then using numerical methods construct 
the additional world volume flux which restores it. 

To check the supersymmetry condition (3.5), the pull-back of the NS-NS two form field 
B 2 = b(p)f)i can be obtained by calculating 

Vi = yi\(wi=ij,) = -^33 dz3 A dz 3 + 5 3 4 dz 3 A dz± + -B43 dz 4 A dz 3 + B 4 i dz^ A dz± , (4.12) 



10 



where the various components are: 

#33 = M {Z1Z3-Z3Z4) (l+ -rW) + -^(^1 + ^2) - ^r( z i ~iz 2 ) } , (4.1:-!) 



/ z 3 z 4 \ ■■>,,■> : ■"•' " 



B 34 = M O4Z3 - 2 3 z 4 ) ( — - \ - i(z x z 2 - ZxZi) + (\z 3 \ + \z 4 \ ) 7^~( z i + iz i) 7^( Zl _ iz2 ^ 



B 43 = i \ (z 4 z 3 - z 3 z 4 ) I -j— pj- I + *(«i^a - Z1Z2) - {\z 3 \ 2 + \z 4 \ 2 ) + —£-(zi + iz 2 ) + — J-(«i - ^2) 



(4.14) 

- 2 )}, 

(4.15) 

B 44 = i I (z 4 z 3 - z 3 z 4 ) f 1 + -j — |g- J JtT^ 1 + iE ^ + ~7^P^ Zl ~ iz2 ^ I ' ( 4 - 16 ) 

Notice that 771 is now only invariant under the U(l) subgroup, and the invariance can be 
made manifest by changing into the coordinates {103, 104,103,104}. 

The components for the pull-back of the Kahler form J = /C'(p)r/5 + K."(p)f}± can also be 
obtained by calculating: 

V5 = V5\(vji=h) = ^33 ^ z 3 A dz 3 + J' 3i dz 3 A du + J' i3 dz± A dz 3 + J 44 dz 4 A dz 4 , (4.17) 
?74 = ?74|(ui 1 = / u) = ^33 dzz A d^3 + J 34 dz 3 A d^4 + J 4 g dz± A ^23 + J^ dz^ A dz4 . (4.18) 
Here the various functions are given by: 

jt , , N^ t/ _-, , NP w _f3|4 , _£4£3 ,,-qx 

"^33 ~~ 1 ' |,,|2 ' ^44 — x "i" I..I2 ' ^34 ~~ |,,|2 ' ^43 ~~ |,,|2 ' V^- 1 ^ 

IMI I Ml I mi imi 



03 - ^3 ^777|2 > 54-^4 ^— 12 • (4.21, 



^33 ~~ S3S3 , J34 — S3S4 , J 4 g — 5453 , J 44 - — S4S4 , (4.20) 

(e 2 - (zl + zp)z 3 __ (e 2 - (fg + zj))z A 

~2\l4 2 ' 2^ 

Similarly, 774 and 7)5 are now only invariant under the U{\) isometry rotating the phase 
of W3/W4. Given the expressions for B2, and the pull-back of the Kahler form J for the 
embedding (4.10), we can now calculate the wedge products explicitly and check if the 
supersymmetric condition (3.3) is satisfied. Let us calculate them in turn: 

- a ~ .(23^4-24^3) /, ,2 , I |2 , I |2 , I |2\ jr. 

r jb f\i-i 1 =i — l\zi\ +\z 2 \ + \z 3 \ +\z4\)dil 

Imp 

|e| 4 sinhpcoshpsin (^±M sin (^=M 

= j-^ '- ^ J -dn (4.22) 

ImT 

with d£l = dz 3 A dz± A dz 3 A dz A , and we have used the explicit coordinates on the deformed 

conifold given in Appendix A. Notice that the overall expression in (4.22) is in fact real as 

i(z 3 Z4 — ZiZ 3 ) = 2Im(z4Z3). Similarly 

i(z 3 u ~ Z4Z 3 ) 



m A 7)1 = 4 ^ («i^- - ZjZi) 2 dO, 

|e| 6 sinh 3 psin (&&a\ sin f^=^ 



dO . (4.23) 

ImI 2 
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Here we notice that the summation in (4.23) is in fact positive definite, since each term 
in the summation (including the negative sign) is nothing but (Im(ziZj)) 2 . Combining 
(4.22) and (4.23), we can finally write down the expression for J AB2 (including the scalar 
functions): 

3 A B 2 = ^) H °i 'H ^ J sinh2 , (4±^M) * , (4.24, 

3|/i| z \smh(2/3j — 2p/ 

where the scalar functions KL'{p) and b{p) are as given in (2.32) and (2.31). The expres- 
sion (4.24) is manifestly non-vanishing and we conclude that in the absence of additional 
worldvolume flux F2, supersymmetry is broken. 



4.3 Restoring Supersymmetry with Worldvolume Flux 

Because B2 is not primitive, it is clear that D7-branes in the u>-embedding are not su- 
persymmetric without the addition of worldvolume flux. In this section we will attempt 
to construct the necessary flux. After setting up the problem, we explore the asymptotic 
behavior of the flux in sections 4.3.1-4.3.4 before turning to a numerical solution in section 
4.3.5. The results are summarized in section 4.3.6. Part of our solution agrees in appropri- 
ate asymptotic regimes with a previous proposal in [35], with some corrections which we 
describe in detail. 

Our strategy is to first find a basis of 2-forms on the D7 worldvolume which are (1,1) 
and primitive. Then we take a linear combination of these forms and impose the Bianchi 
identity, which gives a system of partial differential equations for the coefficients of the 
basis forms, and we proceed to solve this system. 

To make the calculations more tractable we will suppose the parameter /j, is much 
greater than the deformation parameter e so that we may work in the KT region of the 
deformed conifold. In this region, the background geometry is simply that of the warped 
conifold with no deformation but with nontrivial background fluxes. The advantage of 
taking this limit lies in the symmetries of the background. In taking w\ = \i on the 
deformed conifold the isometries of T 1,1 are broken down to a single U(l), but in the KT 
limit the isometry is U(l) X U(l). In addition to the U{\) on the deformed conifold which 
rotates w% and W4 by opposite phases, there is an extra U{\) in the KT region which acts 
as W2 — ► e 2 * Q iL>2,u>3 — > e ia W3,W4 — > e JQ u>4. The U(l) 2 isometry implies that instead of 
having a four-variable problem, we can use the isometry to eliminate the azimuthal angles, 
reducing our problem to a problem of two spatial variables. 

It is useful to express our embedding in terms of explicit coordinates on the conifold. 

fsinf 

2/3 



Upon restricting to W\ = r 3 ' 2 e 1 ' 2 (^ 4>i M gj n ?x s [ n ?z = ^ we see th a t we m ay eliminate 



^37 , V< = <Ai + 02, (4.25) 



and take p to be real. In terms of these angular coordinates, large r corresponds to Oi 
and when 0\ = 62 = vr we are at the minimal radius to which the D7-brane extends. 
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When pulled back to the locus w\ = fj,, the ttij in (2.6) may be written as 

„ o dvjA A dwA 

On = -2z/i 2 * 2 4 4.26 

n Q . 2 dw 3 Adw 3 

r. , r. o • 2 w 3' iJ] id'WA A ^3 + WA_W 3 dw 3 A (^4 

iZl2+ii 2 l = -22// -; n tj 13- 2W1 I2~~i 2T~ (4.28) 

|W3| |W4|(|W3| + I 1 ){\ W M + A 4 ) 

To check primitivity, we need the Kahler form (up to a factor) pulled back to the D7 brane 
worldvolume : 



J oc Qifiu + Q 2 22 + cot y cot y (O12 + 2 i) (4.29) 



where 



Q i = ^ + cot 2 ^. (4.30) 

There are two simple anti-self-dual (1,1) forms which are relevant to us: 

Xi = fin + 1 tan ^ tan ^ Q 2 (Jfci + ^12) (4.31) 

X 2 = fi 22 + - tan -± tan yQi (0 2 i + ^12) (4.32) 

There is also a third anti-self-dual (1,1) form, but it plays no role in the later calculations: 

X 3 = d6i A d6 2 - sin B\ sin 6> 2 #i A d<fc. (4.33) 

It is not hard to see that any component of X 3 in T will either induce a violation of the 
Bianchi identity or will be singular. 

We can take the following linear combinations of X\ and X2 which have nice properties: 

P = QxXx - Q 2 X 2 = Q1O11 - Q 2 ^22 , (4.34) 

Q = X 1 -X 2 . (4.35) 

The form P is closed, dP = 0, while Q has a particularly simple exterior derivative: 

lift rift 

dQ = -F^- A 22 - — 4" A On • ( 4 - 36 ) 

sm Pi sm #2 

We will use P and Q as basis two-forms for the worldvolume flux, and search for an 
anti-self-dual form JF = a(Oi,02)P + P(0i,02)Q which satisfies the Bianchi identity. It is 
enlightening to write the forms P, Q explicitly in complex coordinates: 

P (3 \w 3 \ 2 \ dw 3 A dw 3 _ /3 \wa\ 2 \ dw A A dwA 

"•" ..2 n... 12 1 ,.2\2 o + ..2 n„.. 12 1 ,.2\2 V 4 -J/J 



2i/i 2 V 2 V 2 J (\w 3 \ 2 + li 2 ) 2 \2 y? J(\wa\ 2 + h 2 ) 2 
Q dw 3 A dw 3 dw4 A dw. 



2ifi 2 (\w 3 \ 2 + fi 2 ) 2 (\wa\ 2 + ijl 2 

J ~V4 



fi 2 fi 2 \ wsWAdwA A dw 3 + WAiv 3 dw 3 A dw4 



W3P l^4| 2 / (l^3| 2 + ^ 2 ){\Wi\ 2 +/i 2 ) 
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(4.38) 



By counting powers of 1^3,4, we see that in the limit 11)3,104 — » 00, Q is non-singular while P 
is potentially log-singular. On the other hand, when W3 or u>4 vanishes, P is non-singular 
while Q can be log-singular. 

If one prefers to parameterize the D7-brane world volume by the angular coordinates 
(#i, 0i, 02, ^2)) the basis forms may be written as 

'3 



P 



Q 



(71 

- + cot 2 — 
2 2 



d0\ A sm.0\d§\ 



— h cot 2 — ) d6*2 A sin 02U.fi 1 



(4.39) 



1 , 01 , 02 

- tan — cot — 
2\ 2 2 



cot — tan — ) (d02 A sin 0\d<f>\ + <i#i A sin 02t 



1) 



(4.40) 



+d6>i A sin 0]d(/>i - d^2 A sin 2 <#2- 
Let us now turn to the Bianchi identity which T must satisfy. The B-field in the 



Klebanov-Tseytlin geometry is Bi = Sgs 2 a log ^1^2, so that pulling back H3 = dB2 to the 



D7-brane world volume we have 
dp = H 3 = - 



9sMa> ( Ji 2 

— - — I cot —d0\ + cot — d0 2 



AW2 



(4.41) 



Plugging in the ansatz T = a(0\, 02)P+(3(0i, 02)Q, we obtain a first-order system of partial 
differential equations for <x{0\,02) and /3(#i,#2) : 



de^a = -— d fl2 /? - yr-de^ + . (3 + -=- cot -± 
Q2 Q2 Q2SW.01 Q 2 2 

de 2 a= yz-de.P - —d d2 P + . — /3 + ppcot^. 
Qi Qi Qism^2 Qi 2 



where we have 
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COS 02 — COS 01 



g s Ma' 



(4.42) 



(4.43) 



sin 0i sin 02 4 

This system of differential equations is challenging to solve but we will show that 
numerical methods combined with some analytic tricks will allow us to find a solution. 
The system with two functions can be converted to a second order partial differential 
equation for the function (5(0\,02) by differentiating and eliminating dg 1 dg 2 a. Explicitly, 
the second order equation of interest takes the form 



de 1 -—do 1 + 



Q 



Q2 1 



h h 



+ 



+- 



sin #2 



1 
Ox 



1 



1 

Q~2 

1 

Q~2 



h + 



1 



Qi sin 2 

(3(01,02) 



de 1 



1 



Q2 sin 0\ 



-k ( cot — dg 1 



cot 



-o, 



'«■> 



1 

Q2 



(4.44) 



smc/i 

1 

.QiV " 2 

Our main task in the following subsections will be to solve this complicated equation. 

Once we have solved for /3(9i,02), we can then find the second function a(#i,#2) by 
the following procedure. From the system (4.42), we construct the quantity (5| + dg 2 )a. 
This procedure results in a standard Poisson-type equation for a(0i, 02) with a complicated 
source term which depends on the function (5(0\, 02)- 
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4.3.1 Boundary Conditions and Finiteness 

Before turning to the numerical analysis of the flux equations, we need to identify the 
boundary conditions which we expect the function /3(0i,#2) to satisfy (Recall that the 
angular coordinates are defined over the range (0, ir).) We claim that there is a unique set 
of physical boundary conditions on the supersymmetry-restoring flux. For the function /3, 
the boundary conditions turn out to be compatible in a certain large-radius limit with the 
proposal of Benini [35] . 

The primary physical constraint on (3 is that we expect its contribution to the flux, 
(3Q, to be finite. One then finds that the allowed boundary behavior of (3 is tied to the 
asymptotic properties of the two- form Q. There are three important regimes to consider: 

• 9± = tc or 62 = vr: When one or the other of the angles 9i = tt, the only physical 
choice of boundary conditions is Dirichlet: 

/?(7T,0 2 ) =/?((?!, 7T) = 0. (4.45) 

The simplest argument for this comes from examining the basis two- form Q in (4.38). 
For example, in the limit 9\ — > tt, we have {w^l — > so the basis form Q exhibits a 
logarithmic divergence. To exclude this singular flux, we must set (3 = 0. 

• 02 — > with 9\ > 0, or vice versa: In this limit, \w/^\ — > oo with \w%\ finite, and 
one sees from (4.38) that Q is 0(1). Thus the only constraint in this limit is that (3 
is finite. 

• $1)02 — ¥ simultaneously: In this limit, Q actually vanishes as l/|w| 2 , with 
1^3,4 1 —* oo. Then finiteness of the flux only requires that (3 grows no faster than 

. \ a as 9i — > 0. 

sin 0i 

To specify the large radius (small 9) boundary conditions precisely, we will now study the 
second and third limiting cases in detail. The second case may be usefully regarded as a 
"factorization" limit in which the D7-brane world volume splits into two branches. Along 
one branch we have 9\ — > and along the other, we have #2 - ► 0. The two branches connect 
when both #i,#2 — > 0, corresponding to the third case, which represents an "interpolating" 
limit between the two branches. 

For the second function a, the associated asymptotics may be studied similarly. The 
basis form P is regular in the interior of the D7-brane worldvolume, but has logarithmic 
singularities when 9\ or 92 vanishes. This logarithm is physically acceptable, as it is 
compatible with the log growth of the fluxes in the KT solution. Therefore when either 
angle vanishes, we require that a grows no faster than logr. On the other hand, in the 
interior of the brane, we expect a to be finite. 

4.3.2 Large Radius "Factorization" Limit 

Let us take the limit #2 - > 0, with 9\ held fixed. The various terms in (4.44) have the 
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behavior 



S-tan|l (4.46) 

Q2 - ^ (4.47) 



Multiplying the equation for (3 by 02 and taking the limit as 62 — > 0, we see that (4.44) 
reduces to the following, discarding terms of order Q\: 



d 9l — tan -±d 9l + - tan -j9 2 d0 2 9 2 d e2 - d 9l Q-feSft, + ^-d dl + d 01 ( — J 



li 



= -2kd 6l (i-) (4.48) 

Unfortunately this differential equation is not separable, but it does have homogeneous 
scaling with 62 - all terms are of zeroth order. 

In the previous subsection we argued that in this limit (5 should be finite. Thus the 
leading small 62 behavior should be for (5 to be constant as a function of Q 2 . The equation 
(4.48) simplifies to 

ddi ^ tan y 5ei/3 + dei (qt) = ~ 2k (4 - 49) 

and the solution is 

= -2k + °\ ( cos O1-8 log sin ^ ) + ° 2 (4.50) 

2sin 2 |- V 2 7 2sin 2 ^ 

with two integration constants which we need to determine. We may fix one constant 
by requiring consistency with the Dirichlet boundary condition (3{9\ = ir) = 0, which is 
satisfied if 

-2*- | + | = °- ( 4 - 51 ) 

We will be able to obtain a second constraint after we have performed the analysis in the 
next subsection. 

Now we can substitute into the first order equations to solve for a. In the small 62 
limit these reduce to 

d 9l a = -^tan ^-e 2 de 2 P (4.52) 

tan^f 1 1 2k 

e 2 de 2 a= - 7 ^d dl [3-—e2d e2 P+—(3+—. (4.53) 

Vl VI Vl VI 

The first equation reduces in our limit to dg-^a = because (5 is independent of 62. The 
second equation with (5 plugged in gives a ~ — 2c\ log #2- 
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4.3.3 Large Radius "Interpolating" Limit 

When 61 and 62 are simultaneously small, the asymptotic analysis of section 4.3.2 breaks 
down. For example, there can be terms of the form ^ which we discarded in the strict 62 — > 
limit but are order one when the two angles are small simultaneously. This immediately 
raises the question of whether it is possible for the solutions along the two branches 0\ = 
and 62 = to be mutually consistent when both angles are small, so that the flux smoothly 
interpolates between the two branches. In the following analysis we answer this question 
in the affirmative. 

To proceed, we have found the following change of variables to be useful: 

/ 9 8 \ 
u = — log I sin — sin — I (4.54) 

(sin ^ \ 
s4J (455) 

where u is clearly related to the original radial coordinate r as in (4.25) The natural limit 
to take with both of the Oi small is simply 

u — > 00 , v held fixed. (4.56) 

Now we turn to the decoupled second order equation (4.44) for /?, in the u and v 
coordinates. We may expand the resulting equation in powers of e~ u for large u, in the 
form 

(Co + e~ u O l + e~ 2u 2 + ...)/3 = g + e~ u 9l + ... (4.57) 

with homogeneous and inhomogeneous parts separated. It is tedious but straightforward to 
find (multiplying the equation obtained from (4.44) by an overall factor to reduce clutter) 

Oq = -rr^+2 cothw— + 1 (4.58) 

ov z ov 

cosh-u / d 2 d d 2 \ 3 , d 2 

Ox = — : — ( 37^-2—- tt^ ) +-sinhu- — 



4 \ du 2 du dv 2 J 2 dudv 

— (7sinhuH — J- 2coshv (4.59) 

2 \ smhvj ov 

go = ~2k (4.60) 

g\ = 4fecoshv (4.61) 

In principle one can now expand (3 in powers of e~ u and construct a solution iteratively, 
order by order in e~ u . 

At the zeroth order in the expansion, we encounter a small surprise: in the large u 
limit, the leading homogeneous part of the equation for (3 reduces to an ordinary differential 
equation, 

OoP h = 0. (4.62) 
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This equation is easily solved: 

Ph = /i(«)^- + /2(«)^4- (4-63) 

smn u smn w 

Demanding that (3 is non-singular at v = sets /2 = 0. 

At leading order, fi(u) can be any function of u, but in the full solution its form 
must be compatible with all the other asymptotics we have identified. This eliminates 
most possibilities for f\. For example, a natural ansatz one might try is to set (3 = 
—2k + /i(u)iEEv "^ 0( e ~ u ) with fi(u) of zeroth order in e~ u , so that the leading homoge- 
neous solution and the leading particular solution are of the same order. However, a little 
manipulation shows that this ansatz must be wrong, because the homogeneous part van- 
ishes when v — ► ±oo. The resulting function cannot possibly satisfy the Dirichlet boundary 
condition (4.45). 

What we need is a form for the function f\{u) so that the homogeneous piece fth 
approaches a constant in the limit u — ► oo, v — ► ±oo (which corresponds to (6\, 62) — ► (0, ir) 
or (7r, 0). There is only one such function which suffices, up to subleading terms: 

Ph = be u (-^—)+... (4.64) 

\ u smn v J 

Note that the growth of this (3h for u — > 00 is sufficiently slow to be physically acceptable, as 
determined in section 4.3.1. At the next order in e~ u , we find that the system of equations 
is solved by 

= p h - 2k + 0(1/m) (4.65) 

We will fix the coefficient b in the next subsection. 

Let us now turn our attention to the second function a. The coupled first order 
equations for a and j3 in the u and v variables, expanded to leading order in e~ u , are: 

8 u a = -e~ u (cosh v {d u (3 + (3 + 2k) + 2 sinh vd v /3) (4.66) 

d v a = -e~ u (cosh vd v (3 + sinh v{(3 + 2k)) (4.67) 

To leading order at large u (and including the first subleading term in 1/u), this is solved 
by 

a = a - 2blogu - b VC ° th - + 0(l/u 2 ) (4.68) 

u 

Because the basis form P was closed by itself, a can always be shifted by a constant ao- 

One curious feature of the function (3 which is clear from the large-it expansion is 
that the first subleading term for small #2 depends as l/log(sin -j); in radial coordinates 
this appears as l/log(r) rather than as a power of r. The function l/log(sin -y) has the 
amusing property that although it vanishes when 62 — > 0, its first derivative with respect 
to 62 diverges. We will make use of this subleading behavior in the next section. 
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4.3.4 Consistent Boundary Conditions and Subleading Asymptotics 

Having found asymptotic solutions in two different large radius regimes, we now must check 
that these solutions are mutually consistent. Along the way, we will need to compute the 
next-to- leading asymptotic behavior of the solutions to our PDEs. 

The first additional consistency requirement will fix the ambiguity we left in (4.50). 
There we had two integration constants, c\ and C2, with one relation given by demanding 
consistency with the "small radius" Dirichlet boundary condition. To fix the second inte- 
gration constant, we can take the simultaneous limit Q\ — > and #2 — > (recall that we 
have solved for the exact 6\ dependence in the small 62 limit) and then compare with the 
results of the "interpolating" asymptotic; the results should agree. For the limit of small 
02, we can take (4.50) with (4.51) imposed to find 

^-(2 t + cl )(l--^)-J^logsin| (4.09, 

whereas for (4.65) written in the angular coordinates, we find 

26 / log sin ^\ 

^^li 1 - 2 ^]- 2 ' (4 ' 70) 

(In (4.69) and (4.70) the ~ is meant to imply that subleading terms of order ( log sin -y I have 

1 ■ e i 
been dropped, although we have retained a term proportional to $-.) Demanding that 

log sin -y- 
lo sin ?!■ 

the terms proportional to a constant times — —g- and — g 2 are in agreement sets 1 

sin -^ sin -J* 

d = (4.71) 

b = k. (4.72) 



Note that with c\ = 0, we have 



/? = 2k cot 2 — (4.73) 



in the leading small 62 limit. This asymptotic for (3 are identical to the proposal of [35]. 
The interested reader will find a comparison of our notations in Appendix D. 

To obtain a more detailed matching, and to see the leading nontrivial behavior of a, we 
must go on to the first subleading order. Expanding in the natural expansion parameter, 
l/log#2) we write 



01 | h 

log 02 log 2 e 2 



P=-M + [3o + T £h r + TT^TT + ■ • • ( 4 - 74 ) 



1 One slightly subtle point is that in the factorization limit 82 — > 0, in principle we should also collect 
contributions at all orders of e"'°^"', from the point of view of the interpolating asymptotic. However, the 
additional contributions which we drop are not proportional to — — j— so they do not affect the matching 

sin 2 -J- 

of these particular terms. 
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The function (5q was given in (4.73) and the solution for (3\ has the same form: 

A = 2i%(™ Sl - 8l08Si 4) + 2 iW- <« 5 > 

As with the leading order, we set the two integration constants as follows. One constant 

is set by demanding consistency with the small radius Dirichlet condition and the second 

i ■ e i 
is set by comparing with the term proportional to g- — 2 e in (4.70) (the other terms 

log sin ~y sin 2 -± 

cannot be matched directly without working at higher order in the large-u expansion, 
because we have discarded 1/u corrections in (3h which can shift these terms.) These two 
constraints set 

C3 = c 4 = k. (4.76) 



so 

/3i = ^(cos 2 |-41ogsin|) (4.77) 

Now, to find the leading behavior of a, we plug into the the equations (4.52), (4.53) 
and solve. The answer is analytic but somewhat complicated: 



a = ao — 2/clog 



1 ^ 
log y 



k ( . o . Ox n . 0i , 0! 

- — 27T ~ lo § sin — + 2 log sin — log cos — 

log &2 ^ A ^ *■ 

+ 2logsin| + ^Li 2 (sin 2 |) 



The leading log log is consistent with the behavior in the large-it limit of (4.68), but its 
appearance is perhaps a bit surprising from the perspective of the "factorization" limit 
in Section 4.3.2. In the factorization limit, at leading order the log log is not necessary to 
solve the field equations (although it could have been included) , but it is actually necessary 
to include it to solve consistently at first subleading order. 

4.3.5 Numerical Solution 

In this subsection we construct a numerical solution for the functions a and (3 characterizing 
the D7-brane worldvolume flux; the numerical computation gives strong evidence that the 
D7-branes in the w-embedding are supersymmetric with the addition of this flux. 

In the preceding analysis, we identified a self-consistent set of asymptotic solutions 
for the differential equation (4.44). The task that remains is for us to show that these 
asymptotics correspond to a smooth solution over the entire worldvolume of the D7-brane. 

Under more fortunate circumstances, we would be able to invoke existence and unique- 
ness theorems for solutions of this PDE, and the task of showing that the D7-brane is 
supersymmetric would be complete. As it happens, the mathematical problem we must 
solve is not of a standard type. The second-order PDE for (3 is elliptic through most of its 
domain, but it fails to be elliptic on the line 9\ = 62, where the equation vanishes identically 
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(recall that the second order equation was obtained by differentiating and subtracting the 
first order equations.) We are not aware of general theorems for this type of equation 2 . 

In lieu of an abstract existence proof, we will attempt to demonstrate the existence of 
a solution numerically, employing a finite element method. Strictly speaking, this method 
is only guaranteed to converge for elliptic problems, but we will proceed anyway under the 
assumption that near-ellipticity is enough. The danger that we face is that along the line 
9\ = 02 where ellipticity fails, everything is a solution, including spurious behavior, and 
the numerics can potentially be quite unstable. 

We have used the publicly available FreeFEM++ package [36] to perform the numerical 
calculations. This software uses the variational, or "weak" formulation of the finite-element 
problem. We use the default UMFPACK algorithm. In the calculations shown below, we 
have used a mesh with 2558 vertices, and we have checked by varying the mesh that the 
solutions we obtain are not mesh-dependent. 

A trick that is useful in taming the spurious solutions of our non-elliptic PDE is to 
impose Dirichlet boundary conditions whenever possible. Fortunately the asymptotic anal- 
ysis of the preceding subsections allows us to do precisely this. The boundary conditions 
we have chosen are that j3 vanishes when either angle 9i = n as in (4.45) and along the 
boundary where Q{ — > we require that (3 is given by the asymptotic form in (4.65) with 
b = k as determined by the analysis in section 4.3.4. 

Although the flux due to [3 is finite, in our formulation the function j3 diverges as 0i 
and 62 simultaneously vanish, so we have solved instead for a rescaled function 7 defined 
by 

2 01 2 02 



7(0i,0 2 ) = I 1 -cos 2 ycos 2 - 1,3(01,02). (4.79) 

which does not have this singular behavior. We also impose a small cutoff, 0j > 9 C = 10 -7 
to prevent problems with division by zero. To be specific, our boundary conditions for 7 
take the form 

7 (0 1 ,7r) = 7(7r,0 2 ) = O (4.80) 

.„ „ , 2k log sin % - log sin % 

7 (0i, C ) = -2k+ . . I , . I (4-81) 

siir ^ log sin ^ + log sin -^ 

,„ „ s , 2k log sin %- log sin % 

7(0 C ,0 2 ) = -2k+ . . I / . I (4.82) 

sm z ^ log sin ^ + log sin ^ 

This form was chosen to agree with the asymptotic solution (4.73) in the limit C — > and 
to be compatible with the solution (4.65) in the "interpolating" limit. 

The result for 7 is shown in Figure 1. We see that the gradient of 7 appears to grow 
near the boundaries 9% — > 0, consistent with the subleading asymptotic dependence 1/ log 0j. 



2 Incidentally, the fact the we obtained an ordinary differential equation rather than a PDE in the limit 
of section 4.3.3 was a sign of this degeneracy. Sometimes it is possible for such a degenerate equation to be 
rendered elliptic by cleverly dividing by zero, but in this case it is not clear to us that such a technique is 
applicable, as it would make the asymptotic ODE (4.62) singular. 
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Figure 1: The 7 component of the flux, as found by a numerical calculation. We have set k = 1. 

Near the origin, this effect becomes more pronounced; we believe it reflects the behavior 
of the subleading asymptotics in the Oi coordinates and does not represent a physical 
singularity. The most important thing to take from this diagram is that the solution exists 
and it appears to be free of spurious behavior. 

For the second function a we have implemented the procedure of constructing a Poisson 
equation with a /3-dependent source term. The boundary conditions on a are of generalized 
Neumann type - the normal derivatives of a must be compatible with the first order system 
of equations (4.42). Strictly speaking, these boundary conditions leave our problem slightly 
ill-posed, but the ambiguity simply corresponds to shifting a by a constant, which was 
expected because the basis form P was closed. Figure (2) shows our result for a. The 
mildly singular behavior at the origin appears to be compatible with the log u dependence 
expected from (4.68), and near the corner (#1 = it, 62 = it) it appears that a flattens out 
and approaches a constant. 

4.3.6 Summary of Results 

Because the calculation of the supersymmetry-restoring flux was quite complicated, let us 
pause to summarize the result. The self-dual (1,1) flux can be written as a linear combi- 
nation of two basis two-forms P and Q which are defined in (4.35). We have constructed 
a full solution numerically which appears to be smooth everywhere in the interior of the 
brane worldvolume, and which is presented graphically in Figures (1) and (2). 

The fact that our solution for the D7-brane flux is numerical is somewhat unfortunate, 
as it makes the flux awkward to use in calculations. However, the numerical solution 
does have relatively simple behavior at large and small radius, and we have given analytic 
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Figure 2: The a component of the flux, computed numerically in terms of the angles #j. We have 
set k — 1. 



expressions for the large radius asymptotics. The leading asymptotics are consistent with 
the proposal of Benini [35] ; we have extended his result by demonstrating that the flux can 
be continued over the entire D7-brane worldvolume. 

The small radius asymptotics {0\ ~ n and 02 ~ tt simultaneously) are that (3 = and 
a is a constant. 

When one angle (say 62) is taken to be small, the leading behavior of (3 is 



/? = 2£;cot 



2 p i 



(4.83) 



which is consistent with the earlier proposal of [35] (see Appendix D), and a has the leading 
behavior 



a = ao — 2/clog 



log 



(4.84) 



Here the log log behavior gives a deviation from the result of [35] . 

When both angles are small, the asymptotic analysis is somewhat different. In terms 
of the coordinates defined by (4.54) and (4.55), the leading asymptotic is given by 

1 v 



(3 = ke u 



u sinh v 



2k 



(4.85) 



where we have included the —2k at the first subleading order in e u , and the leading 
asymptotic for a is 



a = ao — 2klogu. 



(4.86) 
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4.4 Other Embeddings 

In this subsection we point out that several well-known examples of D7-branes embed- 
dings are not supersymmetric without the inclusion of worldvolume gauge flux. It would be 
interesting, though presumably quite difficult, to study whether a supersymmetry-restoring 
flux exists in each case. 

Interestingly, given a holomorphic embedding there is a simple criterion that in many 
cases can immediately show that D7 worldvolume flux is necessary for supersymmetry. 
When the D7-brane is embedded holomorphically, it is possible to write its embedding 
equation in terms of the Wi coordinates: g{wi) = 0, say. Then if g{wi) is a polynomial, in 
the limit of large radius r one may approximate g{wi) by truncating to its terms of highest 
rank. If the resulting homogeneous polynomial contains a factor of one of the Wi (up to 
an SO(4) rotation), then in the large radius limit the brane embedding admits a branch 
which locally looks like the simple u>-embedding, which as we have already seen requires a 
nontrivial flux to be supersymmetric. 

Let us now turn to some examples that have been studied in the literature. The 
holomorphic embedding (4.1) admits a generalization that was proposed in [27] 

za = I (zf + zl) , (4.87) 

where f(x) is an arbitrary function of x = z\-\-z\. However, we may also write x = 2(w\W2), 
so if f(x) is a polynomial, then in the large radius limit / ~ {w\W2) p for some p, and it 
may be approximated by p copies of the w\ = and W2 = embeddings. Based on the 
argument given above, we would expect this embedding to be non-supersymmetric without 
a worldvolume flux, and we have checked this explicitly in Appendix B. A similar argument 
applies to the Karch-Katz embedding [25]: 

2 2 

w lW2 = Z A±^ =/ / } p G C , (4.88) 

and the interested reader will find details of the check in Appendix C. Finally, there is an 
interesting class of embeddings proposed by Arean et al in [37] given by 

4 

J]>f=^ P ( 4 -89) 

This embedding admits a factorization in terms of the w% in the large radius limit, so based 
on the examples we have studied so far we expect that it will require a worldvolume flux 
to be supersymmetric. 

One might hope that, as in the case of the w-embedding, that the necessary worldvol- 
ume fluxes might be constructible for more general embeddings - the fact that the flux we 
constructed vanishes for the part of the D 7- brane that dips into the throat is a promising 
sign. The explicit calculations are daunting, however. A general holomorphic embedding 
will break all the isometries of the conifold, and the resulting system of equations one would 
need to solve to find the flux can be a system of partial differential equations for three func- 
tions of four variables. Ideally one would hope for a more abstract proof of existence for 
solutions, and we leave this as an interesting but difficult open problem. 
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5. Field Theory Remarks 

The D7-branes in the warped conifold geometry that we have studied in the previous 
section have an interesting interpretation in the dual gauge theory. In this section we make 
some remarks about the field theories for the various brane embeddings we have studied 
in the earlier part of the paper (For a good review of cascading field theories, see [38].). 

When no D7-branes are present, the dual gauge theory of the warped deformed conifold 
is a non-conformal field theory with a product gauge group SU{N+M) x SU{N) and matter 
fields A\ t 2,Bi t 2 which transform in the bifundamental color representations (N + M, N) c 
and (N + M, N) c . The theory also has a superpotential 

W = ATr ■(A i B j A k B l )e ik <? 1 . (5.1) 

When M = 0, the gauge theory is actually superconformal [39]. Moreover, one can show 
that the moduli space of the field theory corresponds to the moduli space of N D3 branes 
probing the conifold geometry. To see this, take some number of D3-branes and move 
them away from the tip of the conifold by giving expectation values to the fields Ai,Bj. 
Imposing the F-term and D-term equations, one discovers that with the identifications 

wi = AiBi, w 2 = A 2 B 2 , W3 = AiB 2 , u>4 = A 2 Bi . (5.2) 

we recover precisely the defining equation of the conifold (2.7). 

When M/0, conformal invariance is broken and the theory undergoes renormalization 
group flow [19, 31]. Along the RG flow, the two gauge group factors take turns being 
strongly coupled; when one of them goes to strong coupling, one can perform a Seiberg 
duality which results in a weakly coupled description. Each Seiberg duality causes the gauge 
group to change from SU(N + M) x SU(N) to SU(N) x SU(N - M) [40]. Ultimately, in 
the far IR, this "cascade" must end, as the gauge group ranks cannot be negative. At that 
point, the theory undergoes confinement, and has a dynamically generated scale which is 
related to the deformation parameter of the conifold [19] . 

When we add D7-branes in the gravity theory, the gauge theory is correspondingly 
modified by the addition of a number of matter fields charged as fundamentals of the 
gauge group. The introduction of additional matter should change the duality cascade 
pattern, and it was argued in [28] that the extra matter causes the duality cascade to slow 
down as the RG scale decreases. 

5.1 Cascade Pattern for z embedding and its Generalizations 

The analogue of the KS duality cascade when we include probe branes embedded by 
z\ = /j,, is actually quite simple, and is manifestly self-similar, as was pointed out in [41]. 
Although the discussion in this subsection will contain nothing new for experts, for the 
sake of explicitness, we review here the cascade analysis. 

Because the worldvolume flux can vanish for the supersymmetric z-embedding, it is 
reasonable to guess that the appropriate field theory description for the z-embedding in the 
KS background differs from that of the singular conifold theory only by changing the ranks 
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of the gauge groups from SU(N) x SU(N) to SU(N + M) x SU{N). The superpotential 
for the case of one such D7 brane (which one can motivate based on RG flow from a related 
orbifold theory [28]) is 



W Z1=I1 = X iq {A a B a - fi)q + \ 2 ^j^ . (5.3) 



the additional quarks q and q transform in the (N + M, 1) and (N + M, 1) representations 
of the SU(N + M) x SU(N) gauge group (of course we can also choose q and q to be 
charged under the SU(N) gauge group instead.) One check of this superpotential is that 
if we probe the theory with a D3-brane, then the quarks, which one can think of as D3-D7 
strings, become massless precisely when the probe D3-brane intersects with the D7. The 
superpotential has an obvious generalization to the case of K coincident D7-branes. Let us 
consider the case where the additional quark q transforms in the (N + M, 1) representation, 
so that the quark is coupled to the SU(N + M) factor which flows to strong coupling. As in 
the case of the regular KS cascade, there are mesonic operators M a p = A a Bp, and because 
of the added quarks there are also operators which we label as C, D and N: 

{C a )™={A a )™if a , (D a ) a m = q?(B a y m , K = q?qt, i = l...N+M, a,b=l,...,K. 

(5.4) 
These additional SC/^iV + M^-invariant quantities and along with the original mesons M a p 
allow us to rewrite the total modified superpotential for the z embedding into: 

TV 2 
W Total = hTr (M a pM lS ) e a ^ 5 + XiTr(D a C a - fiN) + A 2 — . (5.5) 

Now, as the SU(N + M) factor flows to strong coupling, we see that we have N c = N -\- M 
and Nf = 2N + K, so we are in the conformal window and may perform a Seiberg duality. 
Under this duality, the weakly coupled theory with gauge group SU(N—(M—K)) xSU(N), 
and 2N + K flavors coupled to the SU{N — (M — K)) factor. Notice that in contrast with 
the original KS theory, where the difference between the ranks of two gauge groups remain 
unchanged throughout the duality cascade, here the difference between the ranks reduces 
by K at each step of the cascade. The resulting dual superpotential now receives the 
following additional contributions: 

- \Tr(M a ^a a b^) + Tt{C a q'b a ) + Tr(D a a a q') + Tr{Nq'q')} . (5.6) 

Here a a and b a are the dual fields for A a and B a , under the SU(N — (M — K)) x SU{N) 
gauge group, they transform as (N — (M — K), N) and (N — (M — K), N) representations 
respectively; whereas the K q' and q's, the dual quarks and anti-quarks for q and q trans- 
forms in the (N — (M — K), 1) and (N — (M — K), 1) representations. 

Now, we may integrate out the SU(N — (M — K)) singlets. First, note that one obtains 
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the following F-term equations: 

M Q/3 : 2/ie a7 e /3<5 M 7<5 + -a a \P = , (5.7) 

? 

N : \ 2 N - /J, + -q'q' = , (5.8) 

? 

C a : AiA* + \q'b a = , (5.9) 

D a : AiCc^-a^ = 0. (5.10) 

Substituting (5.7)-(5.10) into (5.5) and (5.6), we obtain the superpotential for the weakly- 
coupled SU(N -{M-K))x SU{N) theory: 



W DuaI = -£2 ^lV(o Q 6^a 7 6 5 )e^e^ + ^q 1 (b a a a - ^) q' + ^(q'q') 



(5.11) 



This superpotential (5.11) has exactly the same functional form (up to field re-definitions 
and extra constants) as the original superpotential (5.1) and (5.3), except that the quarks 
are charged under the gauge group which flows to weak coupling. 

Now let us note that if the additional quarks q, q are coupled to the gauge group that 
flows to weak coupling, the analysis proceeds essentially as before, except that it is simpler 
because the quarks q and q do not bind to form mesons. After performing a Seiberg 
duality on the strongly coupled gauge group, the two gauge groups switch roles again and 
the quarks q, q will again be coupled to the gauge group which flows to strong coupling. 
Thus after two steps of the cascade, the theory is fully self- similar. 

5.2 The vacua for Kuperstein-like Embeddings 

As was shown in [27] and reviewed in Section 4.1, the simple z-embedding is supersym- 
metric for all \x on the full deformed conifold. Therefore it is interesting to study the case 
of a D7-brane embedded by z\ = [i in the extreme infrared (/i ~ e), where the conifold 
becomes deformed. In this case, the quarks due to the D7-brane have masses of the order 
of the confinement scale and can modify the IR dynamics of the theory. 

Suppose the number of units of 3-form flux, after having decreased from the D7 back- 
reaction, is M* 3 . To investigate the moduli space of the theory, we use the standard trick 
of probing the theory by adding a single mobile D3-brane. The gauge group on the D3 
brane probe is SU{M* + 1) x U(l), and we have matter fields a a and b a transforming 
respectively in (M* + l,l) and (1,M* + 1) representations, moreover we also have an 
additional quark q' coming from 3-7 string transform in (M* + 1,1). By self-similarity, 
the UV superpotential on the probe here is essentially identical to (5.11), but with the 
fields and the parameters modified accordingly. At low energy, the UV superpotential can 



3 The explicit value M* can be deduced from the pattern of the flavor cascade described earlier to be 
M* = y/{M - K/2) 2 - INK - K/2. 
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be written again in terms of the SU{M* + 1) invariant quantities and it also acquires a 
non-perturbative Affleck-Dine-Seiberg superpotential [42] : 



^A D s = (M*-2)^_ j . (5.12) 

Here the determinant of the 3x3 meson matrix detO is given by 

detft = —f^e&MapMys - e a ^e ps C a D p M l5 , (5.13) 

whereas the various gauge invariant quantities are defined similarly as in (5.4). From (5.12) 
and (5.13) we can deduce the .F-term equations: 

M ap : e^e ps {hM aP - S(Q)(NM a p - C a D p )) = , (5.14) 

N : X 2 N + -S(n)e a7 e^M al3 M y5 = Ai/i , (5.15) 

C a : Xt^Dp = S(n)e^e^ 5 D p M 7S , (5.16) 

D a : \ie a(3 C a = S(n)e ai e p5 C a M l5 ■ (5.17) 

with 



S(U)= li—fr] -j^r. (5.18) 

v ; IdetfM detft v ; 

One class of simple solutions to these F-term equations is given by setting C a = D a = 
and we have 

1 



N =2x~\ X ^ L ± V ( A ^) 2 " 4A 2/idet M af3 



I A 3M* \ M*-2 

hdetM a0 = — — — - - . (5.19) 

p yNdetMapJ y ' 

The deformed conifold branch naturally appears when we consider the limit (Ai/i) 2 <C 
4A2/idet M a ^, which yields: 



A 2 A 6M * 



i 

2M*-1 



detM^^I^^I . (5.20) 

Another extreme limit which can be taken is that (Ai//) 2 3> 4A2^idet M a p, and again simple 
algebra shows that we recover the deformed conifold. These limits precisely correspond to 
moving the probe D3-brane away from the flavor D7 branes, and locally it detects the 
deformed conifold geometry. 

It would be very interesting to understand the other regions of the moduli space especially 
with C a ,Dp ^0. 
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5.3 ^-embedding and its Generalizations 

Now let us move to the cascade analysis for w-embedding and its generalizations, such 
as the ACR embeddings [37]. We shall begin by relating the superpotential for various 
ACR-embeddings to the one for the simplest w-embedding (4.10); we then on the cascade 
pattern for the simplest w-embedding, as the more complicated cases can be analyzed 
similarly. 

Recall that for the w-embedding (4.10), an appropriate additional term to the KS su- 
perpotential (5.1) was given in [28]: 

W vn=ll = Xq(A 1 B 1 -fi)q, (5.21) 

where q can transforms in the representation (N + M, 1) and q transforms in the (N + M, 1) 
representation of the gauge group. This superpotential can be generalized to many other 
interesting cases. Take for example the Karch-Katz embedding [25], w±w 2 = fi- In the large 
radius limit, we may approximate this embedding equation by wxw 2 = 0, which factorizes 
into two branches w\ = and w 2 = 0. On the field theory side, this naturally corresponds 
to introducing two sets of additional quarks {q\} and {q 2 } with the superpotential 

Kigi(Ai.Bi)gi + K 2 q 2 {A 2 B 2 )q 2 ■ (5.22) 

At smaller radius, however, the two branches meet and at least one set of quarks will be 
massive. We can describe this with superpotential terms of the form 

7i<?2<?i + 72<7i<72 with 71 > 72 . (5.23) 

Taking the larger of the two masses and integrating out the associated quarks, we obtain 
a superpotential of the form q[(AiBi)(A 2 B 2 ) + 72] 9 ; as desired. Of course if we continue 
further and integrate out quarks using the 72 mass term we will be left with a trivial 
superpotential. 

Now, for the more general ACR embedding, we can implement a generalization of the 
procedure we have described for the Karch-Katz embedding. Again, the ACR equation 
can be factorized into different branches in the vanishing \i limit, each of which is given by 
Wi = for some i, motivating superpotential terms of the form 

Pi P2 PS Pi 

^ K n l qn 1 (AlB 1 )q ni +'^2 Kn 2 Qn 2 (A 2 B 2 )q n2 + /j K>n 3 gn 3 (A\ B 2 )q n3 + ^ K n4 q n4 (A 2 B 1 )q rl4 , 

7ll=l 712 = 1 ^3=1 714 = 1 

(5.24) 
which are simply linear combinations of the simplest tu-embedding. By adding appropriate 
mass terms and integrating out quarks as they become massive, it is not hard to see that 
we can flow to a superpotential of the form 

W ACR = Xq [(A 1 Biy i (A 2 B 2 y 2 (A 1 B 2 )P s (A 2 B 1 )P i - fi p ] q, (5.25) 

which is related to the ACR embedding equation in a direct way. 
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Now, in the case with no three-form flux on the conifold, the D7 branes in the w- 
embedding factorizes further into two branches, using the relation (5.2) w\ = A\B\. A 
superpotential reflecting this factorization property was proposed in [28] and studied in 
[35]: 

W flav = giQA x q + g 2 qB 1 Q . (5.26) 

By introducing the quark mass terms miqq + m 2 QQ , m 2 3> m\ the familiar superpotential 
term (5.21) (up to appropriate redefinition) can be recovered by introducing the quark mass 
terms m\qq+m 2 QQ and integrating out appropriately. The cascade pattern for this type of 
superpotential has been studied in [35] , with the result that the self-similarity under Seiberg 
duality transformation no longer presents in this case and a sequence of additional terms 
involving singlet states and quarks is generated. Explicitly the resulting superpotential 
after a step of cascade is given by: 

AiTr(a Q ^a 7 6 5 )e Q7 e /3<5 + \ 2 {Pa 2 q' + qb 2 P + Z qq) . (5.27) 

Here Ai,2 are combinations of coupling constants and dynamical scale; So is an additional 
gauge singlet field. In the proposal of [35], the cascade is almost self-similar; at each 
step the term with the singlet E acquires additional A, B fields and (at least naively) 
becomes more and more irrelevant. This unusual cascade deserves further study; perhaps 
the worldvolume flux we constructed in section 4.3 may shed some light on the nature of 
the field theory. 

6. Discussion 

We conclude with a few comments on extensions and possible applications of our results. 
Having seen that worldvolume fluxes can restore supersymmetry, it would of course be 
interesting to study fluxes for embeddings other than the simple linear u;-embedding, and 
to extend our calculations to the deformed conifold. As we have already remarked, this is 
in general a daunting mathematical exercise, but perhaps it will be possible to prove an 
existence theorem of some kind for the necessary fluxes, possibly based on the topology 
of the embedding. Understanding the worldvolume flux may also help in extending the 
calculations of D7-brane backreaction along the lines of [41, 43]. It would also be interesting 
to study D7-branes in warped cones other than the conifold. 

In the gauge/gravity duality, these D7-branes add flavor to the dual KS gauge the- 
ory, with corresponding prospects for understanding mesons and baryons in a confining 
theory. It might be interesting to revisit the calculations of the meson spectrum, as for 
example in [44], with worldvolume fluxes turned on. In particular, it would be interesting 
to understand better the physical meaning of the closed (1,1) form P which we identified 
in Section 4.3, which one can turn on without breaking any supersymmetry; it should 
correspond to deformation of the field theory by some operator. A clearer understanding 
of these supersymmetry-restoring fluxes may also be important in clarifying the nature of 
the duality cascade with probe w-D7-branes. 

There are also potential applications to inflationary model building. The model con- 
structed in [22] used a very special brane embedding to generate the inflationary potential, 
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but with several different supersymmetric brane embeddings, one can imagine combining 
D7-branes in various ways to obtain many other inflationary models. It would also be very 
interesting to understand if the presence of extra flux can modify the explicit expression for 
the gaugino condensate given in [21]. It was argued in [22] that at the leading order expan- 
sion of D7 DBI action, the additional supersymmetric flux should not modify the results 
obtained in [21] . Combining the severe constraints of holomorphy of the superpotential and 
the global symmetries preserved by the embedding, one is tempted to conclude that the 
potential due to the gaugino condensate in the deformed conifold should be quite similar 
to that of the singular conifold. Having a supersymmetric D7 configuration including the 
necessary world volume flux should allow us check such a statement explicitly, following 
the calculations in [21]. 

Also, it would also be interesting to consider small deviations from the extremal limit, 
and allow supersymmetry to be softly broken by a D-term potential whose magnitude is set 
by J s J A J- where £4 is the four cycle which D7 brane wraps on in the deformed conifold 
[45] (see also [46] for an interesting construction in the resolved conifold.) Combined with 
the fact that £4 is generally warped, this may allow us to consider a variety of interesting 
scenarios for supersymmetry breaking. 
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A. Deformed Conifold Coordinates 

The embedding coordinates {z{\ transform as 4 under the SO{A) isometry group. An 
alternative parametrization for the deformed conifold constraint (2.12) can also be given 
by a 2 x 2 complex matrix: 

w= lU + iZi Z!~iZ 2 \fw 3 W 2 \ 

V2 \z\ + iz2 -Z3 + 1Z4 I \ w\ Wi J 

defining equation for the conifold can then be expressed in terms of W as 

e 2 
det W = W3W4 — w\W2 = — — • (A. 2) 

The coordinates {wi} now transform as (2, 2) under SO(A) ~ SU(2)l x SU(2)r isometry 
group. The complex embedding coordinates of deformed conifold {z\, z 2 , 23, Z4} can be 
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expressed in terms of the real coordinates {p G M , tp G [0, 47r] , #1,2 G [0, 7r] , </>i,2 G [0,27r]}, 
S = /) + iijj as: 



zi = e 



22 = e 



23 = e 



2:4 = e 



cosh ^-j cos ^^^J cos [-^) +ismh {- ) cos [~^ ) ** ^^^ 
-cosh ^-j cos ^^^J sin ^^^J +iBmh ^-j cos ^^^J cos ^—^ 

, /s\ . f fli + M ( <h-<h \ , ■ . ^M ■ / %-M . ^1-^2 

ooahl-Jsmf— g— Joob(— 2— J+.amhf-JBnif— ^— jsmf— 2— 



, ,'5\ . (e x + e 2 \ . (<h-fa 



. , , ,'£\ . / gl~g2 \ / 01 -fe 

' S11 2 / Sm I 2 / COS I 2 



(A.3) 



.3/2 



In the large radius limit p — > 00, ecosh(p/2) ss esinh(/o/2) « £7=-, the deformed conifold 
coordinates (A.3) and (A.l) smoothly reduce to those of singular conifold. 



B. Generalized z-Embeddings 

We next examine a generalization of the holomorphic embedding (4.1) proposed in [27] 

z 4 = f {z\ + zl) , (B.l) 

where f(x) is an arbitrary function of x = z\ + z\ = 2(u;i«;2). This class of embedding 
generally breaks the SO (4) isometry group down to the U(l) which rotates w\ and w 2 by 
opposite phases. 

The general expression for the pull-back of NS-NS B-field in this case is given by: 

171 = r/4| 24= j( x ) = Budzi A dz\ + B 12 dz\ A dz 2 + B 2 \dz 2 A dz\ + B 22 dz 2 A dz 2 , (B.2) 

with the various components given by the expression in (4.2): 

B xl = 2 {(z 3 z 2 - z 2 z 3 )(f'zi - fzi) + (f 2 zi - f 2 zi){z A z 2 - z 2 z 4 )} , 

B 22 = -2 {(z 3 zi - ziz 3 )(f'z 2 - fz 2 ) + (f 2 z 2 - ] 2 z 2 )(z±z\ - Z\Z/$\ , 

B 12 = (z A z 3 - z 3 z A ) + 2 {(zi(z 3 zi - ziz 3 )f + z 2 (z 3 z 2 - z 2 z 3 )f) - (z 3 -> z 4 , f -> f 2 )} , 

B 2\ = («3*4 - Z 4 Z 3 ) - 2 {{zi{z 3 Zi - Z\Z 3 )J + Z 2 (z 3 Z 2 - Z 2 Z 3 )f) - (z 3 -> Z 4 , /' -> / 2 )} . 

(B.3) 



Here 



Z3 = ^ 2 -(^ + / 2 (^)), /'(*) 



df(x) 



f2 



l + f'(x) 



(B.4) 



cte 2z 3 

the expression for the previous case (4.1) is recovered by setting f'{x) = 0. The pull-back 
of the Kahler two-form for the generalized Kuperstein embedding can be written as: 

7)5 = ?75|(^ 4 =/(a;)) = J'i\ dz\ A dz\ + J' 12 - dz\ A d^2 + Jgj dz 2 A dJ?i + J^ dz 2 A d^2 , (B.5) 
7)4 = ^4|(z4=/(x)) = ^11 dz\ A dzi + J^2 dz\ A dz2 + Jgj dz 2 A d^i + J^ dz 2 A d^2 . (B.6) 
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with the components 

J' n = 1 + \zi\ 2 p(x) , J22 = 1 + \z 2 \ 2 p(x) , J' 12 - = ziz 2 p(x) , J' 2l = z 2 zip{x) , (B.7) 
Jil = S±Si , J 2 2 = S2S2 , J12 = 5*1 5*2 , J 21 = 5^2 "^l j (B-8) 

P (x)=4|/(x)| 2 + |1+ | 2/ . / 2 (x)|2 , (B.9) 

1*3 1 

51 = - [(zgf! - Zl z 3 ) + 2z 1 (z 3 za - z 3 )f(x)] (B.10) 

Z3 

52 = — [(z 3 z 2 - z 2 z 3 ) + 2z 2 (z 3 Z4 - z 3 )f'(x)] (B.ll) 

23 

Notice here that these expressions reduces to the ones for simple Kuperstein embedding 
(4.1) when f'(x) = 0. To check the SUSY condition for the generalized embedding (4.87), 
we can first calculate 

fc A VI = {Zl ~ Z2 ~, Z 2 2 ~ Zl \ [z 3 (l + 4) f{x) - z 3 (l + zl)f'(x)] 

1*3 1 

+ (z 3 z 4 - Z4^ 3 )[2|/ / (x)| 2 (l + \z 3 \ 2 ) + f{x) + f(x)])dn , (B.12) 

where dfl = dz\ A dz2 A dzi A dz~2- For an arbitrary f(x), the above expression does not 
vanish unless f'(x) = 0. One can also in principle calculate 7)4 A 771, and demonstrate that 
it is non- vanishing for arbitrary /(x). In fact, one can consider a special case f{x) = x, 
i.e. Z4 = z 2 + z\\ for each specific value of 24, this is identical to Karch-Katz embedding 
(4.88), which we will also show momentarily to be non-supersymmetric. We thus conclude 
that the embedding (4.87) cannot be supersymmetric in the deformed conifold for arbitrary 
function f{x). However, with additional world volume flux, the supersymmetric condition 
can still in principle be satisfied. 

Let us remark also on a statement in [27] that the K-symmetry conditions must be 
satisfied because of a symmetry argument. The argument used the following logic. Under 
the interchange Z\ — > .22,-22 - ► z±, the Kahler form J and the embedding equation are 
invariant, while B2 gets a minus sign. Thus JAB acquires a minus sign under the inter- 
change. There is also a rotational symmetry (from the SO{3) preserved by the embedding 
equation) rotating z\ and 22, under which J, B, and the embedding equation are all in- 
variant. One can write J A B = (j>(zi, z±,Z2, z~2)dz\ A dz2 A dz\ A dz~2- Then (ft is invariant 
under the rotation but acquires a minus sign on the interchange, and [27] claimed that the 
only function which could satisfy both properties is (ft = 0. This claim is not correct, as 
z\Z2 — 22-21 is a counterexample. 

C. Karch-Katz Embedding 

We consider next the holomorphic embedding given in [25] 

2 2 

w 1 w2 = ^^ = fi 2 , f.eC, (C.l) 

which explicitly breaks the SO (A) isometry group into SO (2) x SO(2), and these SO (2) 
subgroups act by rotating the phases of the ratios w\/w2 and w 3 /w4. In the asymptotic 
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limit r 3 S> |/i| 2 limit, the embedding (4.88) can in fact be viewed as two copies of the D7 
brane w-embedding we have earlier (4.10) w\ = and W2 = 0, and they intersect and fuse 
together at some finite radius by turning on D7-D7 interactions. 

The pull-back of the NS-NS two form field in this case is given by calculating: 

Vi = Vi\( Wl w2=!J. 2 ) = B il dz\ A dz\ + B X3 dz\ A dz 3 + B 3 \ dz 3 A dz\ + B 33 dz% A dz 3 (C.2) 
where the various components are given by: 

B x i = \z 2 \~ 2 (z 1 z 2 - z 2 z 1 )(z 3 z 4 - z 3 z 4 ) , 

B lg = - <^ (z 2 z 4 - z 4 z 2 ) + (z 4 z 4 - z 4 zi)— + (z 2 z 3 - z 3 z 2 )^- + (ziza - z 3 z~i)-^- } , (C.3) 
L z 2 z 4 z 2 z 4 ) 

s 3i = \ ( Z 2Z 4 - z 4 z 2 ) + (z 4 z 4 - z 4 z~i)~ + (z 2 z 3 - z 3 z 2 )— + (z x z 3 ~ z 3 z 4 )^-^ \ , (C.4) 
I z 2 z 4 z 2 z 4 ) 

s 33 = \z 4 \~ 2 {ziz 2 - z 2 zi)(z 3 z 4 - Z 3 Z 4 ) , (C.5) 

with 

z 2 = 2fi 2 - z\ , z\ = (e 2 - 2/j 2 ) - z\ . (C.6) 

The components for the pull-back of the Kahler form J can also be calculated from: 

V5 = Vb\(w 1W2=y ?) = J'n dzi A dz\ + J' 33 dz 3 A dz 3 , (C.7) 

7)4 = vA{w-LW2=fj?) = J'il dzi A dz\ + J" 3 dz\ A dz 3 + J 3 ^ dz-$ A dz\ + J 33 dzs A dz% .(C.8) 
The various functions here are: 

J d = 1 + S> 43 = 1 + S' ^3 = 4 = o. (c.9) 

j„_ _ {Z\Z2 ~ ZlZ 2 ) 2 „ _ (Z3Z4 ~ Z 3 Zj) 2 

11 I 19 ' 33 I 19 ' 

\Z2\ 2 \Zi\ Z 

T „ {z\Z2 ~ ZiZ 2 )(z 3 Z4 1 -Z3Z4) „ (Z\Z2- ZIZ2)(Z3Z 4 - Z 3 Z 4 ) ,„ , 

The SUSY condition (3.5) can be checked by calculating the wedge products in turns: 

- a - {ZIZ 2 - Z2Zl){z3Z4~ Z A Z 3 ) , 2 2 2 2 , 

15A)J1 = ; r^ r^ ( \Zl\ +\Z 2 \ + \z 3 \ + \z 4 \ ) dil , 

\Z2\ \Za\ 

where dQ = dz\ f\dz 3 l\dz\ Adz 3 . In (C.ll) we have a product of positive definite sum with 
a generically non-vanishing number — 4Im(zi Z2)^-to.{z 3 z/ 4 ) . Similarly we calculate 7)4 A f]\ 
and obtain: 

[Z\Z2 - Z 2 Z 1 ){z 3 Z A - ZiZ 3 ) Y,i+j=l( Z i Z 3 ~ Z jZi) 2 ,^ 
77 4 A r/l = 1 pr- — p^— — dil . 

\Z2\ \Z4\ 

Once again we have a product of a positive sum and a non-vanishing number. We can 
again show that the SUSY condition (3.3) is not satisfied for Karch-Katz embedding (4.88) 
without additional magnetic field F2. This result is perhaps not so surprising given the fact 
that Karch-Katz embedding can be regarded as two copies of tw-embeddings asymptotically, 
and we have shown that u)-embedding cannot be supersymmetric unless additional world 
volume flux is turned on. 
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D. Comparison with Benini's Proposal 

Benini [35] has proposed a form for the asymptotic flux in the limit of large radius. The 
purpose of this Appendix is to compare our notation with his so that it is clear how our 
proposal compares with his in the relevant limit. 

In the large radius limit, we may take one of the angles 8\ or 82 to vanish. A helpful 
piece of intuition is that the worldvolume of the D7-brane in this limit "splits" into two 
branches, corresponding to one or the other 8i vanishing. 

Let us take the limit where one of the angles, say 82, is small, and then work order-by- 
order in 82- The D7-brane worldvolume in the compact directions can be described by the 
coordinates 9\, 82, (f>i, 4>2- However, to compare with the large radius limit we will trade 82 
and (p2 for r and i/j. 

In this limit, the (1,1) basis forms P and Q take the form 



P 



3 ^ „ dr c 

-fin + 3— A g- 
2 r 



8 dr 

+ cot—(d8i A3 5 -3— AsinMcM 
2 r 



1 8 dr 

Q tan— (d0i As 5 -3— Asin0id0i) (D.l) 

2 2 r 

where g 5 = dip + cos8id(f)i + cos #2^2- The claim in [35] was that the term in square 
brackets in P was a solution to the inhomogeneous Bianchi identity in the large radius 
limit. This corresponds to a combination of P and Q with a = Ik and (5 = 2fccot 2 =&•, or 
in other words 

F = 2kP+ l^ 7r -2k)Q (D.2) 

V s m 2 1 J 

which we have written in a form to make clear the comparison with formulas such as (4.69). 

If one took the large radius limit on the other branch, with 8\ small, the form of T is 
almost the same, with 8\ and <j)\ replaced by 82 and </>2 everywhere, and multiplication by 
an overall factor of — 1 (due to the antisymmetry of the basis two- form 102 which appears in 
the background 2-form potential P>2.) In our formalism, this sign flip is built into the basis 
forms P and Q which are antisymmetric under the interchange of indices 1 <-> 2, while the 
functions a and /3 are symmetric. 

Finally, note that our asymptotic form for a in the large radius limit (4.84) in the 
radial coordinates is 

a = —2k log log r + const + 0(1/ log r) (D-3) 

with a leading log log term; this term is effectively constant in the leading 1/logr approx- 
imation, but it does affect the field equations at subleading order (and therefore has to be 
included), as we argued in Section 4.3.4. 
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